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ABSTRACT 


Theoretical  and  empirical  information  necessary  for  the  design  of 
an  electrical  propulsion  test  facility  collector  system  complex  is  com¬ 
piled.  Design  considerations,  which  include  sputtering,  beam  accom¬ 
modation,  secondary  emission,  and  less  serious  secondary  effects,  are 
considered  in  some  detail. 

A  theoretical  model  is  presented  which  accounts  for  the  energy 
dependence  of  the  sputtering  phenomena  first  at  low  ion  energies  and 
secondly  in  the  region  where  the  Bohr  screening  parameter  is  of  the 
order  of  unity.  Based  on  these  functions  obtained  from  the  theoretical 
model  in  the  two  regions,  an  asymptotic  empirical  relation  for  the  sput¬ 
tering  ratio  as  a  function  of  ion  energy  is  found  to  give  a  satisfactory 
account  of  available  data.  In  addition,  equations  are  developed  to  evalu¬ 
ate  the  effects  of  sputtered  particles  upon  vacuum  chamber  pressure, 
the  amount  of  flux  leaving  the  target  that  returns  to  the  engine,  and  the 
time  required  for  the  target  to  become  saturated  with  beam  particles. 

Two  different  collectors,  the  solid  metallic  type  and  a  liquid  metal  sys¬ 
tem,  are  evaluated.  The  advantages  and  limitations  of  each  are  discussed. 
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1.0  INTRODUCTION 


To  establish  a  ground  test  program  for  long-term  testing  of  elec¬ 
trical  propulsion  systems,  one  of  the  major  problems  to  be  resolved  is 
that  of  adequate  ion  beam  collection.  Although  not  a  problem  in  flight 
testing,  adequate  beam  collection  is  essential  in  the  ground  test  facility 
complex  to  ensure  accurate  simulation  of  space  (Refs.  1  and  2). 

An  efficient  collector  system  is  required  to  accommodate  any  new 
material  and  its  associated  energy  introduced  into  the  chamber  from 
the  exhaust  of  an  operating  electrical  propulsion  system.  Means  must 
also  be  furnished  to  isolate  the  engine  proper  from  any  external  phenom¬ 
ena  foreign  to  its  actual  space  environment  such  as  sputtering  effects 
caused  by  beam-target  interaction  and  secondary  electron  emission.  The 
importance  of  these  depends  dn  (1)  the  duration  of  the  test,  (2)  the  ion  cur¬ 
rent  density,  (3)  the  ion  energy,  (4)  the  ion-target  material  combination, 
and  (5)  the  degree  of  space  simulation  required. 

The  purpose  of  this  report  is  to  furnish  information  necessary  for 
the  design  and  evaluation  of  a  collector  system  complex. 


2.0  SUMMARY  OF  SPUTTERING  EXPERIMENTS 


2.1  INTRODUCTION 

When  a  target  material  is  placed  in  the  path  of  a  positive  ion  beam 
which  has  been  accelerated  through  a  small  potential,  no  erosion  effects 
are  observed.  As  the  accelerating  voltage  is  increased,  a  threshold 
value  is  reached  and  target  atoms  become  dislodged  from  the  surface. 

The  threshold  value  depends  on  the  target  and  ion  beam  combination 
(Ref.  3).  This  phenomenon  is  known  as  sputtering.  The  number  of  tar¬ 
get  atoms  knocked  off  or  eroded  per  incident  ion  is  defined  as  a  sputter¬ 
ing  ratio  or  yield.  In  general,  the  sputtering  ratio  increases  with  in¬ 
creasing  ion  energy  up  to  a  maximum  value,  levels  off,  and  then  decreases 
with  a  further  increase  in  ion  energy.  For  low  molecular  weight  ion  beams, 
this  maximum  value  occurs  at  a  relatively  low  energy,  and  a  sharp  de¬ 
crease  in  sputtering  yield  is  noted  immediately  after  maximum  value  is 
obtained  (Refs.  4,  5,  6,  and  7).  With  heavier  molecular  weight  ion  beams, 
the  maximum  value  generally  occurs  at  higher  energies,  remains  constant 
with  energy  over  a  much  longer  range,  and  finally  decreases  slowly  with 
increasing  ion  energy  (Ref.  7).  Figure  1  is  a  schematic  diagram  of  a 
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typical  yield  curve.  Several  variables  are  known  to  affect  the  yield  curve, 
such  as  the  residual  pressure  in  the  test  container,  the  angle  and  energy 
at  which  the  beam  material  strikes  the  target  surface,  and  the  nature 
and  crystalline  structure  of  the  target  for  a  given  beam  material. 


2.2  BACKGROUND  PRESSURE  VERSUS  YIELD 

Several  authors  have  investigated  the  effect  of  residual  gas  pres¬ 
sure  on  the  sputtering  yield  (Refs.  8,  9,  and  10).  It  is  generally  con¬ 
cluded  that  at  relatively  high  pressures  the  sputtering  ratio  is  smaller 
than  at  lower  pressures.  For  a  given  current  density,  the  sputtering 
ratio  increases  with  decreasing  pressure  and  finally  reaches  a  constant 
value  (Fig.  2).  Wehner  has  explained  this  in  terms  of  contamination  of 
the  target  surface  by  the  residual  gas  (Ref.  8).  According  to  Wehner' s 
argument,  this  pressure  dependency  would  be  observed  at  much  lower 
pressures  than  is  shown  in  Fig.  2  if  the  rate  of  surface  contamination 
by  the  residual  chamber  gases  were  comparable  to  the  rate  of  surface 
erosion  by  the  impinging  ions.  The  condition  for  the  sputtering  yield  to 
be  independent  of  the  chamber  pressure  can  therefore  be  written 

4-  V  v  f  << 

Yonts  noted  a  pressure  dependency  at  5  x  10“^  torr  (Ref.  10),  using 
a  current  density  of  only  0.  1  ma/crn^.  Figures  2  and  3  are  experiment¬ 
ally  determined  curves  showing  yield  versus  pressure  and  current  density, 
respectively.  These  observations  coupled  with  Yont's  findings  seem  to 
substantiate  Wehner' s  contamination  hypothesis.  Also,  it  should  be  noted 
that  back  diffusion  of  the  sputtered  particles  plays  an  important  role  in 
reducing  the  number  of  sputtered  particles  when  the  mean  free  path  of  the 
residual  gas  is  much  smaller  than  container  dimensions. 


2.3  CRYSTALLINE  STRUCTURE  OF  TARGET 

The  crystalline  structure  of  the  target  has  been  found  to  have  an 
important  role  in  the  sputtering  process.  Rol  and  his  colleagues  under 
Kistemaker  studied  the  sputtering  rate  of  single  copper  crystals  by  20-kev 
argon  ion  bombardment  (Ref.  11).  At  normal  incidence  and  a  current 
density  of  0.  5  ma/ cm2,  they  found  the  yield  of  6.  4  for  polycrystalline 
copper,  8.2  for  a  (111)  copper  surface,  and  4  atoms/ion  for  a  (100)  cop¬ 
per  surface.  Yields  for  the  single  crystal  samples  as  a  function  of  the 
angle  of  incidence  were  found  to  show  maximum  and  minimum  points. 
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These  peaks  indicate  that  sputtering  yields  have  maximum  values  in  the 
direction  in  which  the  beam  "sees"  a  large  density  of  projected  lattice 
points  and  a  minimum  where  the  crystal  shows  the  highest  transparency. 
Although  the  study  of  single  crystal  sputtering  offers  a  greater  insight 
to  basic  collision  phenomena,  the  present  work  is  concerned  only  with 
polycrystalline  targets. 


2.4  THE  STATE  OF  SPUTTERED  PARTICLES 

In  addition  to  the  rate  of  production  of  free  particles  by  the  sputter¬ 
ing  process,  it  is  desirable  to  know  the  state,  energy,  and  angular 
distribution  of  these  particles  so  that  the  appropriate  steps  may  be  taken 
to  accommodate  all  material  added  to  the  system. 

A  number  of  investigations  have  substantiated  the  fact  that  sputtered 
particles  are  individual  atoms  which  are  largely  neutral,  some  being  in  a 
metastable  state.  One  of  the  more  recent  and  more  comprehensive  of 
these  studies,  made  by  Fogel  (Ref.  12),  found  that  the  positive  ion  sput¬ 
tering  coefficient  rises  with  incident  ion  energy  until  it  reaches  a  narrow 
maximum  in  the  neighborhood  of  20  to  30  kev  and  then  decreases.  For 
the  particular  cases  which  he  considered  (noble  gases  colliding  with 
molybdenum),  the  value  of  the  positive  ion  coefficient  never  rose  above 
0.  01  percent  in  the  energy  range  up  to  about  50  kev.  No  detection  of 
negative  ions  of  the  target  species  was  reported  by  any  of  the  authors 
reviewed.  For  every  case  considered  here,  over  95  percent  of  the  sput¬ 
tered  particles  were  neutral  atoms. 


2.5  EFFECT  OF  IMPACT  ANGLE 

The  effect  of  the  angle  of  impact  of  a  beam  with  the  target  material 
is  another  item  of  importance  in  facility  design.  The  following  two  basic 
aspects  are  to  be  considered  on  oblique  incidence  of  the  beam  with  the 
target: 

1.  The  sputtering  ratio  increases  with  increasing  angle  from 

the  target  normal  for  practically  all  materials  until  a  critical 
value  is  obtained  where  a  decrease  in  intensity  becomes 
noticeable  (Figs.  4a  through  c).  Figures  4d  through  f  indi¬ 
cate  a  similar  trend;  however,  the  experiments  were  not 
carried  out  throughout  the  entire  range  to  indicate  a  critical 
maximum  value,  and  subsequent  decrease,  in  intensity. 

Wehner  observed  for  low  energies  that  materials  exhibit¬ 
ing  low  normal  sputtering  ratios,  such  as  scandium,  titanium, 
and  vanadium,  exhibit  a  marked  angular  dependency,  while 
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materials  exhibiting  a  high  normal  sputtering  ratio,  such  as 
copper,  silver,!  and  gold,  show  a  less  pronounced  effect. 
Almen  and  Bruce  report  the  absence  of  this  distinction  at 
high  energies  (Fig.  4f). 

2.  The  angular  distribution  of  sputtered  particles  is  signifi¬ 
cantly  changed  by  oblique  incidence  of  the  beam  at  low  ion 
energies  only.  At  low  ion  energies  (below  a  few  kev),  the 
distribution  exhibits  a  definite  persistence  in  the  forward 
direction  (Ref.  13)  as  opposed  to  the  symmetric  distribution 
resulting  in  normal  sputtering.  Wehner  and  Rosenberg 
illustrate  this  effect  in  the  distribution  curves  depicted  in 
Fig.  5.  At  ion  energies  in  the  kilovolt  range  or  greater, 
the  angular  distribution  of  sputtered  particles  becomes  less 
dependent  upon  the  angle  of  incidence  of  the  beam  (see  sec¬ 
tion  2.6.1  and  Figs.  6  and  7). 


2.6  ANGULAR  AND  ENERGY  DISTRIBUTION  OF  SPUTTERED  PARTICLES 


2.6.1  Neutral  Particles 


Angular  distribution  of  sputtered  particles  has  been  measured  by 
several  investigators  (Refs.  5,  12,  13,  and  14).  Considering  the  com¬ 
plexity  of  such  a  measurement,  the  correlation  between  investigators  is 
exceptionally  good,  and  the  following  general  conclusions  can  be  made: 

(1)  For  low  energy  normally  incident  particles  (below  3  kev),  an  under 
cosine  distribution  has  been  observed  (Figs.  8a,  b,  and  c).  (2)  For 

intermediate  energies  (approximately  3  to  5  kev),  the  distribution  is 
very  nearly  a  cosine  distribution  for  all  angles  of  incidence  (Fig.  9). 

(3)  For  higher  energies  (above  5  kev),  the  distribution  is  seen  to  become 
more  concentrated  in  the  direction  normal  to  the  target  regardless  of  the 
incident  angle  (Fig.  9). 

Figures  8a,  b,  and  c  give  the  angular  distribution  for  a  number  of  low 
energy  cases  of  interest  in  electrical  propulsion  system  design.  Unfor¬ 
tunately,  no  equivalent  high  energy  data  is  available  for  these  cases. 

Measurements  of  the  energy  of  sputtered  particles  have  also  been 
extensively  reported  (Refs.  9,  11,  15,  16,  17,  and  18),  and  recent  meas¬ 
urements  are  in  favorable  agreement.  Some  of  the  most  elaborate  meas¬ 
urements  of  energy  distributions  for  low  energy  ions  have  been  made  by 
Wehner  where  several  methods  were  used  to  measure  the  complete  energy 
distribution  of  these  particles.  The  following  conclusions  are  evident: 

The  energy  distributions  are  closely  approximated  by  a  Maxwellian  dis¬ 
tribution  for  low  energy  sputtering  and  deviate  slightly  from  a  Maxwellian 
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distribution  at  higher  ion  energies.  The  average  energy  of  sputtered 
neutral  particles  varies  only  slightly  (5  to  10  ev)  with  impingement  energy 
up  to  approximately  3  kev  and  rises  to  a  value  of  30  to  40  ev  for  50-kev 
impingements.  The  average  energy  of  sputtered  ions  is  slightly  higher 
than  that  of  neutrals  (Ref.  29). 

The  velocity  distribution  for  a  number  of  particular  cases  can  be 
seen  in  Fig.  10.  Average  velocity  is  plotted  against  ion  energies  in 
Fig.  11. 

2.6.2  Re-emitted  Beam  Particles 

In  addition  to  the  sputtering  phenomenon,  a  target  can  become  satur¬ 
ated  with  beam  material  and  in  equilibrium  re -emit  beam  material  at  a 
rate  equal  to  that  of  the  beam  impingement  upon  the  target.  Very  little 
work  has  been  conducted  concerning  the  re-emitted  beam  material.  Recent 
studies  concerning  this  phenomenon  indicate  that  a  large  majority  of 
material  re-emitted  from  such  a  target  comprises  neutral  atomic 
particles.  The  more  reliable  measurements  of  the  degree  of  ionization 
of  these  particles  indicate  a  concentration  of  ionized  particles,  both  posi¬ 
tive  and  negative,  of  less  than  one  percent.  Based  on  previous  discussion 
and  information  from  other  sources  (Refs.  19,  20,  21,  and  22),  it  is 
believed  that  if  the  ions  strike  the  target  at  (or  near)  normal  incidence, 
the  energy  of  re -emitted  beam  particles  having  lost  most  of  their  energy 
in  collisions  with  the  lattice  atoms  will  lie  between  that  corresponding  to 
the  localized  target  temperature  and  that  of  the  sputtered  particles.  If 
this  theory  is  correct,  because  of  the  randomizing  collisions  with  lattice 
atoms,  the  angular  distribution  of  such  particles  should  be  in  the  form  of 
a  cosine.  At  oblique  incidences,  however,  a  significant  amount  of  the 
primary  beam  is  re-emitted  with  a  considerable  amount  of  energy,  indi¬ 
cating  very  few  thermalizing  collisions.  This  effect  is  especially  pro¬ 
nounced  when  ions  in  the  primary  beam  are  of  lower  molecular  weight 
than  target  molecules  (Ref.  21). 

In  addition  to  re-emitted  beam  material  and  sputtered  target  material, 
sputtered  contaminant  materials  from  the  surface  of  the  target  have  been 
noted  by  most  investigators.  Both  negative  and  positive  ions  have  been 
observed  and  offer  a  method  for  detecting  what  contaminants  are  on  a 
given  metallic  surface  (Ref.  12).  Coefficients  for  the  production  of  such 
material  are  usually  less  than  3  to  4  percent  and  should  be  of  negligible 
importance,  especially  when  targets  are  under  long-term  bombardment 
under  high  vacuum. 

2.6.3  Reflected  Beam  Material 

A  number  of  investigators  have  conducted  studies  with  the  so-called 
reflected  beam  molecules  (Refs.  23,  24,  25,  and  26).  A  closer  definition 
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of  the  term  is  desirable  to  discern  the  difference  between  the  reflection 
and  re-emission.  Obviously,  any  molecule  which  strikes  the  surface 
must  remain  there  for  at  least  the  finite  time  duration  involved  in  a 
single  collision  with  a  target  molecule  before  it  leaves  the  surface. 
Therefore,  in  the  strictest  sense,  all  molecules  leaving  the  target  sur¬ 
face  are  re-emitted;  however,  reflection  is  a  useful  term  to  apply  to 
those  beam  particles  that  make  only  one  collision  before  leaving  the 
target.  In  this  case,  the  reflected  particle  energy  is  relatively  very 
high. 

The  studies  of  reflection  of  normal  beams  indicate  that  very  few 
particles  are  reflected  for  incident  ion  energies  greater  than  a  few 
hundred  electron  volts.  Most  of  the  studies  were  concerned  with  meta¬ 
stable  atom's  and  ions  because  of  the  ease  of  detection.  The  following 
characteristics  are  generally  observed  for  reflected  particles: 

1.  The  mean  energy  is  relatively  higher  than  that  of  emitted 
materials.  (This  feature  actually  follows  from  definition.) 

2.  The  angular  distribution  is  cosine  in  nature.  This  property 
is  apparent  since  the  de  Broglie  wavelength  of  an  incident 
molecule  is  much  smaller  than  the  irregularities  of  the 
surface  of  any  solid  at  room  temperature,  leading  to  diffuse 
reflection  from  the  surface. 

3.  The  values  for  reflection  coefficients  increase  with  energy 
to  a  maximum  of  3  to  4  percent  at  30  to  40  kev  and  decrease 
with  further  increase  in  energy. 

4.  Normal  reflection  can  occur  only  for  the  case  in  which  the 
striking  molecule  has  a  smaller  mass  than  the  target  molecule. 

Conceivably,  this  phenomenon  could  become  an  important  one  for  the 
collector  system  design  for  beams  of  high  current  density.  The  oblique 
incidence  of  a  beam  increases  reflection  of  ion  beams,  especially  where 
the  ion-target  mass  ratio  is  less  than  unity.  Although  several  investiga¬ 
tors  have  noted  this  effect,  no  qualitative  data  have  been  found  in  the 
literature. 


2.7  SECONDARY  ELECTRON  EMISSION 

In  addition  to  the  ejection  of  atomic  particles  from  a  target  under  ion 
bombardment,  the  phenomenon  of  secondary  electron  emission  has  been 
observed.  There  are  two  basic  types  of  secondary  electron  emission: 

(1)  electron  induced  and  (2)  ion  and  neutral  particle  induced.  The  first  is 
more  conducive  to  experimental  and  theoretical  analysis  and  has  been 
under  extensive  study  for  many  years.  A  number  of  books  and  journals 
cover  this  phenomenon  in  detail  (Refs.  26,  27,  and  28),  and  values  of  the 
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secondary  electron  emission  coefficient  have  been  extensively  tabulated. 
The  importance  of  this  phenomenon  in  collector  system  design  work  will 
depend  upon  the  nature  of  the  beam  electrons  to  be  collected.  Under 
normal  circumstances,  electrons  in  the  beam  oj.  an  ion  engine  are  rela¬ 
tively  low  in  energy  so  that  electron-induced  secondary  electron  emis¬ 
sion  may  be  relatively  unimportant  in  collector  system  design. 

The  second  process,  ion-induced  secondary  electron  emission,  can 
be  expected  to  be  more  important  for  this  application.  Experimental 
data  for  low  energy  ion  beams  in  the  ranges  from  0.  2  to  1  kev  have  been 
given  in  Refs.  28,  29,  30,  31,  and  32  and  for  higher  energies  in  Refs.  16, 
33,  and  34.  Some  of  these  data  are  reproduced  in  Fig.  12  for  cesium 
impinging  on  various  targets.  Here  it  is  shown  that  the  secondary  emis¬ 
sion  coefficient  (a)  increases  from  0.  01  electrons  per  incident  ion  with 
energies  of  a  few  hundred  ev  to  values  of  the  order  of  unity  and  higher 
for  energies  up  to  20  to  30  kev.  Values  for  a  are  plotted  for  other  alkali 
metals  in  Fig.  13.  The  general  trend  seen  in  Figs.  12  and  13  holds  for 
all  cases  of  ion-target  combinations  found  in  the  literature.  It  is  widely 
accepted  that  the  secondary  electron  emission  (a)  is  an  order  of  magni¬ 
tude  higher  than  positive  or  negative  ion  emission  coefficients. 

An  increasing  target  temperature  results  in  a  decrease  in  a.  Most 
investigators  explain  that  this  effect  (shown  for  cesium  on  nickel  in 
Fig.  14)  is  caused  by  the  change  in  surface  condition.  A  cleaner  sur¬ 
face  usually  exhibits  lower  values  of  a.  Surface  condition  has  been 
shown  to  be  an  important  factor  in  all  of  the  emission  coefficients;  for 
example,  it  is  well  known  that  the  work  function  of  an  oxide -coated 
tungsten  surface  is  lower  than  that  of  a  clean  tungsten  surface,  resulting 
in  higher  emission  currents.  The  effects  of  secondary  emission  must 
be  considered  when  any  electrical  measurements  are  taken  at  the  target. 
Since  the  incident  energy  of  the  ion  beams  is  expected  to  be  in  the  kilo¬ 
volt  range,  the  secondary  electron  emission  coefficient  can  be  expected 
to  be  the  order  of  unity  or  greater.  The  necessity  for  electrical  isola¬ 
tion  of  the  collector  from  the  test  article  is  apparent  but  will  not  be  con¬ 
sidered  here. 


2.8  EFFECT  OF  VARIOUS  IONS  ON  YIELD  AT  HIGH  ENERGY 

Alme'n  and  Bruce  (Ref.  9)  used  a  wide  variety  of  ion  beams  to  acquire 
further  insight  into  the  effects  of  sputtering.  Figure  15  contains  the  data 
obtained  by  these  authors  using  70  different  ion  species  impinging  upon 
tantalum,  copper,  and  silver  targets  at  an  incoming  energy  of  approxi¬ 
mately  45  kilovolts.  Sputtering  yield  is  plotted  versus  the  atomic  number 
of  the  ion  beams.  From  this  figure  it  is  seen  that  sputtering  phenomena 
at  this  energy  are  not  simply  mass  dependent  but  exhibit  an  atomic  group 
periodicity. 
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3.0  SUMMARY  OF  SPUTTERING  THEORIES 


Several  theories  have  been  advanced  to  describe  sputtering  phenom¬ 
ena  (Refs.  35  through  41).  These  theories  generally  fall  into  two  major 
classes:  (1)  that  of  macroscopic  evaporation  caused  by  local  thermal 
spikes  and  (2)  microscopic  collision  theories.  The  evaporation  theories 
have  been  discredited  because  of  their  inability  to  account  for  (1)  the 
invariance  of  the  sputtering  ratio  with  target  temperature,  (2)  the  ion 
mass  dependence  of  the  sputtering  ratio,  and  (3)  secondary  emission 
caused  by  ion  bombardment  being  considerably  less  than  would  be  ex¬ 
pected  from  thermionic  emission  at  the  surface. 

An  earlier  victory  for  the  evaporation  theories  stemmed  from  the 
Knudsen  angular  distribution  of  sputtered  particles  that  was  observed  in 
experiments.  Wehner  has  recently  shown  that  the  cosine  distribution 
holds  only  for  ions  of  intermediate  energies.  An  "under"  cosine  distribu¬ 
tion  is  generally  noted  for  low  energy  bombardment  while  an  "over"  cosine 
distribution  (heavier  concentration  is  found  in  the  normal  direction)  is 
witnessed  at  higher  energies  (Ref.  13). 

A  low  and  intermediate  energy  range  sputtering  equation  has  been 
developed  by  Keywell  (Ref.  3  5),  based  upon  a  neutron  cooling  model. 

The  validity  of  the  development  is  difficult  to  ascertain  from  experimental 
data.  However,  the  assumption  of  a  random  walk  phenomenon  does  not 
appear  valid  because  of  the  persistence  of  velocity  effect;  that  is,  random 
walk  theory  assumes  that  there  is  no  preferred  direction  of  motion  after 
collision.  The  same  shortcoming  is  evident  in  other  theoretical  develop¬ 
ments  which  assume  the  applicability  of  the  diffusion  equation  (Ref.  38). 

Rol  (Ref.  3  6)  has  finessed  the  problem  of  predicting  the  number  of 
displaced  particles  that  are  able  to  escape  the  target  surface  by  assuming 
that  the  sputtering  ratio  is  proportional  to  the  number  of  particles  created 
at  the  surface.  In  his  treatment  all  ion  collisions  after  the  first  were 
assumed  negligible.  No  theoretical  basis  exists  for  such  an  assumption. 
Langberg  (Ref.  3  7)  developed  a  low  energy  sputtering  equation  assuming 
a  Morse-type  interaction  potential  between  neighboring  atoms.  He  also 
concluded  that  a  rigid-sphere-type  potential  between  collision  partners 
is  invalid  at  low  energies;  this  conclusion  tends  to  contradict  Bohr's  work 
(Ref.  42). 

Several  authors  have  treated  the  problem  of  high  energy  sputtering  of 
low  molecular  weight  particles  (Refs.  38,  39,  and  40).  These  develop¬ 
ments  are  generally  made  considering  Rutherford-type  collisions  and  are 
concerned  only  with  the  portion  of  the  yield  curve  that  decreases  with 
increasing  energy.  Since  the  Rutherford-type  phenomenon  generally 
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occurs  at  energies  higher  than  is  of  interest  to  electrical  propulsion  col¬ 
lector  system  design,  except  for  very  light  collision  partners,  no  further 
mention  will  be  made  of  the  above  works. 


Bohr  (Ref.  42)  has  extensively  treated  the  penetration  of  atomic 
particles  through  matter.  Since  a  single  model  of  the  penetration  phenom¬ 
ena  covering  the  entire  energy  spectrum  could  not  be  readily  attained,  he 
defined  two  parameters  which  may  be  used  to  ascertain  the  appropriate 
model  that  will  best  approximate  the  collision  mechanism  for  the  energy 
range  of  interest.  These  parameters  are: 


k  =  2  ZiZt 


(la) 


t  _  2  Z{  Zt  Rh  M{  ah _  (lb) 

E  fl  as 

Equation  (la)  is  the  ratio  of  the  collision  diameter  to  the  de  Broglie 
wavelength  of  the  incident  particle;  hence,  its  value  dictates  whether  or 
not  the  wave  properties  of  the  particles  can  be  neglected  and  a  classical 
approach  used.  For  values  of  k!>|>:  1,  the  de  Broglie  wavelength  is  neg¬ 
ligibly  small  and  the  use  of  classical  mechanics  is  completely  justified. 
This  approximation  is  valid  even  for  interactions  involving  low  molecular 
weight  ion-target  combinations,  such  as  helium  and  beryllium,  for  ener¬ 
gies  extending  to  1  mev.  Thus  the  wave  mechanical  approach  is  not 
required  for  a  collector  design  and  will  be  disregarded. 


When  an  ion*  interacts  with  a  target  atom,  the  electron  clouds  of 
each  tend  to  overlap  and  screen  the  fields  of  the  two  nuclei.  To  describe 
the  interaction,  Bohr  has  suggested  use  of  a  potential  function  of  the  fol¬ 
lowing  form: 


u(r) 


61  e  2 
r 


e 


(2a) 


The  factor  as  in  the  above  equation  is  a  screening  constant  given  by 

a  = _ _ 

■(z/3  +  (2b) 

Equation  (lb)  is  the  ratio  of  the  collision  diameter  to  the  above  screening 
constant.  If  the  parameter  of  Eq.  (lb)  is  very  much  less  than  unity,  the 
screening  is  weak,  and  Rutherford-type  scattering  may  be  employed 
(provided  k  >  >1).  If  on  the  other  hand,  £  >>1  the  phenomenon  approaches 
a  rigid  sphere  interaction  in  which  the  scattering  of  the  ions  approaches 
spherical  symmetry  in  the  center  of  mass  reference  frame.  Figure  16 
is  a  schematic  diagram  of  the  logarithms  of  Eqs.  (la)  and  (lb).  The 


^Although  it  is  generally  felt  that  the  impinging  ions  are  neutralized 
before  actual  surface  contact  is  made,  the  bullet  particle  will  be  alluded 
to  as  an  ion. 
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region  in  the  upper  left  hand  quadrant  with  horizontal  hatching  is  that  in 
which  the  screening  is  weak  wherein  classical  mechanics  may  be  used. 
The  region  having  vertical  hatching,  which  includes  the  two  lower  quad¬ 
rants  and  a  portion  of  the  upper  right  hand  quadrant,  is  the  region  which 
must  be  described  in  terms  of  wave  mechanics.  Bohr  has  used  approxi¬ 
mate  methods  to  describe  the  phenomena  in  this  region.  The  open  area 
is  a  region  in  which  exact  wave  mechanic  treatment  must  be  used.  In 
the  region  above  the  45 -deg  line  indicated  by  dots,  the  scattering  of  the 
ion  can  be  treated  classically  and  is  the  rigid  sphere  region,  at  least 
down  to  points  near  the  origin.  In  the  following  analysis  the  sputtering 
phenomenon  is  treated  only  in  this  region. 


4.0  ANALYSIS  OF  SPUTTERING  PHENOMENA 


4.1  DERIVATION  OF  SPUTTERING  EQUATION 

In  view  of  the  uncertainties  and  limitations  inherent  with  existing 
theories,  a  theoretical  model  is  developed  in  the  following  to  evaluate 
the  effects  of  sputtering  on  collection  system  design.  ^Consider  the  dis¬ 
tribution  function  f(p)dp,  which  is  defined  as  the  probability  that  a  pri¬ 
mary  displaced  particle  with  initial  energy  Ep  will  come  to  rest  between 
p  and  p  +  dp.  The  mean  range  of  the  particle  may  then  be  defined  as  the 
following: 

oo 

RP  =  P  =  J  p  f  (p)  dp  (3) 

which  assumes  that  the  range  is  very  much  less  than  the  thickness  of  the 
material  through  which  the  particle  is  traversing.  The  probability 
g(p)  that  a  particle  will  travel  a  distance  p  or  greater  before  coming  to 
rest  may  be  calculated  from  the  above  distribution  function. 

oo 

g(p)  =  /  f  (p  ')  dp'  (4) 

p 

The  probability  that  a  displaced  particle  will  escape  the  surface  start¬ 
ing  al;  a  distance  x  from  the  surface  of  the  plate  pictured  in  Fig.  17  may 
be  written  as  follows: 

ww  -  2 "  /  *  (~^r)  s<« is  (5) 

where  S  (0)  is  defined  as  the  angular  distribution  of  particles  with  a  veloc¬ 
ity  component  in  direction  of  the  surface.  The  2?r  arises  from  the  fact 
that  the  scattering  is  assumed  to  be  symmetric  about  the  incident  direc¬ 
tion.  The  above  distribution  function  will  in  general  be  dependent  upon 
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the  mass  of  the  ion  and  target  materials.  The  total  number  of  primaries 
that  are  able  to  escape  the  target  may  then  be  written  as  follows: 

oo 

Nout  =  /  N(x)  W(x)  dx  (6) 

0 


The  number  of  primary  displaced  atoms  per  unit  path  length  of  the 
incident  ion  N(x)  is  defined  as  the  displacement  cross  section  times  the 
number  of  target  nuclei  per  unit  volume.  For  a  rigid  sphere -type  inter¬ 
action,  the  displacement  cross  section  is  given  by  the  following  expres¬ 
sion: 


=  i t  d. 


y  E(x) 


(7) 


The  quantity  in  brackets  represents  the  probability  of  a  displacing 
event  occurring  per  collision  and  7rd^|  is  the  collision  cross  section. 
Since  the  energy  of  the  incoming  ion  will  vary  with  path  length,  the  num¬ 
ber  of  particles  produced  per  unit  path  length  is  dependent  upon  x,  the 
depth  of  penetration.  Although  E(x)  is  not  generally  known,  the  displace¬ 
ment  cross  section  as  a  function  of  x  may  be  obtained  by  dividing  the 
beam  flux  into  two  separate  groups:  (1)  the  ions  which  have  yet  to  make 
their  first  scattering  collision  and  (2)  the  ions  which  have  made  at  least 
one  collision  with  the  target  lattice.  The  probability  of  a  displacement 
collision  for  the  first  type  particle  is  given  by 

Et 


Pd,  =  1 


Y  Ei 


(8) 


The  probability  of  displacing  an  atom  per  collision  for  the  second  type  of 
ion  may  be  written  as  follows: 


where 


Pd, 


1  - 


y  E 


E  i 

E  =  /  E 


dE 

Ei 


Ei 

2 


(9) 


since  the  probability  of  finding  a  group  2  particle  in  the  energy  range  dE 
is  dE/Ei  for  rigid  sphere  encounters.  The  total  probability  of  a  displace¬ 
ment  collision  per  ion-target  impact  may  be  obtained  by  weighting  each  of 
the  above  probabilities  by  the  fraction  of  ions  in  each  group. 


h~  Ei  1 

e  ~x/\  + 

1  - 

2  Et 

1  -  e“x/X 

_ l 

y  E, 

=  1  +  —  (e"x/X  -  2) 

y  Ei  \  ) 


(10) 
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where  the  exponential  term  is  the  probability  that  an  ion  will  make  a 
collision  between  x  and  x+dx  (Ref.  53).  Thus  N(x)  is  obtained  by  multi¬ 
plying  Eq.  (10)  by  N7rd-jJ.  Equation  (6)  may  now  be  evaluated  using 
Eq.  (10). 

f  f  id  dy  g(y)  S(0)  cos  6 


where 


Nout  =  (277)  77  d12  N 


2  E 


77/ 


y  Ei 


y  Ej 


—  f  f  d0  dy  g(y)  S (6)  cos  6 


77/2  00 

J  f  iOdy 


y  c 


os  0 


g(y)  S(0)  cos  6 


(11) 


cos  0 

The  first  two  terms  of  Eq.  (11)  may  be  evaluated  without  a  knowledge  of 
the  distribution  function  g(y).  For  brevity  let 

tt/2 

Vh  =  2  77  /  S(0)  cos  6  dd  (12) 

Integrating  the  first  two  terms  of  Eq.  (11)  by  parts  (the  y  dependency)  and 
dg 


noting  that  ~  =  -f(y),  we  obtain 


j 


Nout  =  W  dl2-  N  <  V>1 


y  E[  o  o 


J  y  f(y)  dy  - 


77/2 

/  /  d  0  dy 


y  Ei 

y  cos  0 


/  y  f(y)  dy 


g(y)  S(0)  cos  9 


(13) 


The  first  two  integrals  appearing  in  Eq.  (13)  are  simply  Rp  defined  by 
Eq.  (3).  Equation  (13)  can  therefore  be  written  as 


Nout  =  77  d12  N 


■At  Rp  (i 


2Et  \ 
y  e;  J 


+  z  (E;) 


where  Z(Ej)  represents  the  last  term  of  Eq.  (13). 


(14) 


To  evaluate  Z(Ei),  a  knowledge  of  g(y)  is  mandatory.  One  would 
expect  the  distribution  function  defined  by  Eq.  (4)  to  have  a  maximum 
value  at  the  surface  x  =  0  or  y  =  0  since  primaries  produced  at  that  point 
will,  in  general,  be  much  more  energetic  and  hence  will  have  a  maximum 
probability  of  traveling  a  distance  y  or  greater.  For  increasing  values 
of  x,  the  distribution  function  is  expected  to  approach  zero,  since  pri¬ 
mary  displaced  atoms  produced  at  a  considerable  depth  within  the  solid 
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material  have  on  the  average  very  little  energy,  and  hence  their  prob¬ 
ability  of  reaching  the  surface  decreases  with  increasing  depth  of  origin. 
A  distribution  function  having  these  general  characteristics  (y  >  0)  is: 


g  (y)  dy 


(y  +^~) 


~  a  (y  +  -f-) 


dy 


/  (y  +  ur) 


—  a  (y  +  -§-) 


dy 


(15) 


a  n  +  1  / 

=  T  (n  +  1)  (y  + 


a  (  y  + 


— ) 
a  ' 


dy 


The  range  distribution  function  may  be  obtained  by  taking  the  negative 
derivative  of  Eq.  (15).  This  is  seen  to  be 


f  (y)  dy  =  laye 


(y+-f~) 


y  + 


(16) 


where  k  is  a  normalizing  factor. 


Figure  32  compares  the  function  above  with  data  obtained  by  Davies 
and  Sims  (Ref.  43)  for  xenon  on  an  aluminum  target.  The  correlation  is 
satisfactory.  Using  the  above  distribution  function  for  g(y),  Z(Ei)  of 
Eq.  (14)  may  then  be  evaluated  (setting  n  =  2).  Upon  integrating  over  y, 
Z(Ei)  becomes: 


Z  (Ei) 


y  Fh  V  (3) 


n, 4 

f d 0  S  ( 0 )  cos  0 


ab  +  (ab)2  +  (ab)3 


(17) 


where 


b  = 


_A _ 

cos  6 


and 


18  l 
5  Rp 


(18) 


The  relationship  between  a  and  Rp  may  be  obtained  from  Eq.  (16)  and 
evaluating  Eq.  (3). 


The  mean  range  of  particles  through  matter  has  been  theoretically 
evaluated  by  Davies  and  Sims  (Ref.  43),  who  extended  Nielsen's  deriva¬ 
tion  (Ref.  44).  The  expression  they  developed  is 


K  1  .19  Mt2 

(z^3  +  z^)/2 

Mt  +  M, 

L  P  Z i  Zt  J 

(19) 


where  K  is  a  function  of  the  mass  ratio  only.  The  values  of  K  were  calcu¬ 
lated  by  Monte  Carlo  techniques,  and  its  value  for  particles  of  equal  mass 
is  approximately  1.  55  x  10 "6.  For  primary  displaced  particles  moving  in 
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the  target,  Eq.  (19)  reduces  to: 

1.19  K  Mt  (2  Zt)'3Ep; 


.V,; 


Rp  = 


2  Zt  p 


=  K  Er 


(20) 


The  average  energy  of  the  primary  particle  in  Eq.  (20)  may  be  written 
in  terms  of  the  incident  ion  energy.  For  isotropic  scattering  in  the  center 
of  mass  system,  the  relation,  neglecting  E-t,  is 


y  Ei 

Ep  =  /  E 

0 


dE 

yEi 


yEj 

2 


(21) 


Equation  (17)  may  be  written  in  terms  of  the  incident  ion  energy  with  the 
aid  of  Eqs.  (18),  (20),  and  (21)  as  follows: 


Z  (Ei) 


r  o) 


n/2 

f  S  (0)  cos  6 


36 

5 


b 


2  2 
Ky  Ej 


+ 


(-¥-b> 


K  y  Ej 


K3  y  E* 


(22) 

d0 


Each  term  on  the  right  hand  side  of  Eq.  (22)  depends  inversely  on  powers 
of  Ei  and  will  take  part  in  shaping  the  initial  portion  of  the  yield  curve 
(note  the  initial  bend  in  Fig.  1).  These  terms  will  rapidly  approach  zero 
in  the  energy  range  of  interest. 


Since  the  range  of  any  secondary  particles  produced  by  the  primary 
displacements  is  significantly  less  than  that  of  the  primary  particles 
themselves,  it  is  assumed  that  the  number  of  secondary  displaced  atoms 
escaping  the  target  boundary  is  negligible.  Thus,  the  number  of  particles 
escaping  the  target  per  incident  ion,  which  is  defined  as  the  sputtering 
ratio,  is  given  approximately  by  Nout  (Eq.  (14)). 

By  combining  Eqs.  (14)  and  (20),  it  is  possible  to  obtain  an  expres¬ 
sion  for  the  sputtering  ratio  as  a  function  of  incident  ion  energy  (neglect¬ 
ing  the  contribution  made  by  Z(E^): 

SL  =  -rdj  N  i/»!  K  y  [Ej  -  H]  (23) 

where 


Since  Eq.  (23)  was  developed  for  the  case  of  Ei  >i>Et,  the  threshold 
energy  H  corresponds  to  the  linearly  extrapolated  threshold  (Fig.  1). 
Equation  (23)  predicts  a  linear  dependency  of  the  sputtering  ratio  for 
incident  ion  energies  well  above  the  threshold  but  still  within  the  domain 
of  validity  of  rigid  sphere  collisions.  This  linear  dependency  has  been 
substantiated  experimentally  (Ref.  45). 
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As  the  energy  of  the  incoming  ion  beam  is  increased,  Eq.  (lb)  indi- 
cates  that  the  screening  parameter  approaches  unity.  In  this  region, 

Bohr  pointed  out  (Ref.  42)  that  the  major  portion  of  the  deflections  take 
place  in  the  potential  field  of  the  inverse  square  type.  Assuming  a  field 
of  this  nature,  the  collision  cross  section  takes  the  following  form: 

=  JT.|  a°  (24) 

The  displacement  cross  section  may  then  be  obtained  by  multiplying 
Eq.  (24)  byThe  probability  of  a  displacement  collision.  However,  at 
energies  much  greater  than  the  threshold  energy,  the  probability  is 
essentially  unity.  Again,  neglecting  secondary  particles  escaping  the 
boundary,  the  displacement  cross  section,  as  a  function  of  depth  of  pene¬ 
tration  x,  in  the  energy  range  £  ~  1  is  obtained  by  considering  molecules 
that  have  made  at  least  one  collision  separately  from  those  that  have  not. 
The  total  displacement  cross  section  is  then  given  by  weighting  the  cross 
section  for  each  type  event  by  the  fraction  of  particles  within  each  group. 
Substituting  for  £  from  Eq.  (lb)  into  Eq.  (24)  yields 


2  77  e  j  e  t  as  (  M;  +  Mt  )  ah  e  A 
e  Mt  E} 


(25) 


simplifying  Eq.  25 


4  77  ej  et  as  t  Mj  +  Mt  )  Rh  ah 
e  y  Mt  E; 


2  77  ej  et  as  (  Mj  +  Mt  )  Rh  ah 


e  Mt  Ei 


-  1 


(26) 


Thus  the  sputtering  ratio  at  £  -  1  may  be  obtained  from  Eqs.  (26)  and 
(6)  and  the  condition  that  N(x)  =  No-.'  The  maximum  value  of  the  sputtering 
ratio  may  be  readily  obtained  since  S  decreases  with  increasing  ion  energy 
in  the  Rutherford  region.  Again  neglecting  inverse  powers  of  Ei, 

Srnax  =  -2-?-12-  N  ^  [4  77  eietas  (Mi  +  Mt)Rhah1  K  (27) 

IVlt  V  e  L 

where  the  function  ^2  has  been  introduced  to  allow  for  the  possibility  of 
a  different  angular  distribution  function  in  this  energy  region. 


Equation  (2  7)  predicts  that  the  sputtering  ratio  is  independent  of 
energy  within  the  energy  range  where  £  is  of  the  order  of  unity.  From 
the  model  that  has  been  developed,  as  f  varies  from  a  very  large  posi¬ 
tive  quantity  to  the  order  of  unity,  the  asymptotes  of  the  sputtering  ratio 
may  be  obtained  from  Eqs.  (23)  and  (27).  Thus,  the  sputtering  ratio 
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begins  with  essentially  a  linear  rise  beyond  the  initial  bend  and  peaks 
out  at  some  maximum  value  where  £  is  of  the  order  of  unity.  It  has 
been  found  that  throughout  the  transition  region,  the  sputtering  ratio 
may  be  represented  by  the  following  equation: 


S 


=  sr 


-C(Ei-H  ) 
e 


(28) 


where  Smax  is  the  sputtering  ratio  obtained  from  Eq.  (27).  By  expanding 
Eq.  (28)  in  a  Taylor  series,  S  approaches  a  linear  function  of  (E  -  H)  at 
low  energies  as  required  by  Eq.  (23).  Also  at  high  energies,  S  approaches 
a  constant  value,  Smax,  as  predicted  by  Eq.  (28). 


Without  a  prior  knowledge  of  and  ^2,  absolute  values  of  Smax  and 
C  cannot  be  calculated.  However,  to  expedite  the  solution  of  the  collector 
problem,  Smax  and  C  have  been  considered  as  experimental  constants. 
Values  of  these  constants  for  materials  of  interest  are  tabulated  in  Table  1. 


Figures  18  and  19  illustrate  the  fit  of  Eq.  (28)  to  experimental  results 
of  other  workers.  Note  that  the  curves  are  essentially  linear  for  low 
energies  and  approach  a  constant  value  for  £  on  the  order  of  unity,  as  pre¬ 
dicted  by  Eq.  (28). 


4.2  TARGET  SATURATION 

When  an  energetic  ion  beam  impinges  upon  a  target,  sputtering  is 
accompanied  by  other  phenomena  which  result  in  additional  requirements 
for  the  successful  beam  accommodation.  For  example,  the  target  will 
rapidly  become  saturated  with  beam  material  and  subsequently  will  be¬ 
come  a  source  of  free  molecules  of  both  the  target  and  beam  species 
that  rebound  into  the  system.  These  molecules  then  advance  to  other 
points  in  the  system  where  they  are  ultimately  collected.  This  section 
describes  the  effects  of  these  phenomena  upon  test  chamber  conditions. 

Analytical  treatments  to  determine  the  effect  of  terminating  an  ener¬ 
getic  ion  beam  on  test  chamber  conditions  have  been  presented  by  other 
investigators  (Refs.  46  to  48),  but  unfortunately  none  of  these  have  proved 
successful.  Richley  and  Cybulski  (Ref.  -  47)  adopted  a  model  developed  by 
Mickelsen  and  Childs  and  added  two  additional  assumptions:  (1)  they 
assumed  that  the  effects  of  sputtering  could  be  neglected  and  (2)  that  a 
near  unity  collection  efficiency  existed  at  a  target.  *  While  either  assump¬ 
tion  is  valid  in  a  specific  ion  energy  range,  neither  is  generally  valid  in 


^Collection  efficiency  is  defined  as: 

^  No.  particles  leaving  target /  sec 
No.  particles  striking  target /  sec 
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the  range  of  interest.  The  second  assumption  is  valid  only  during  the 
first  few  milliseconds  of  the  operation  of  an  ion  beam  and  is  in  direct 
contradiction  with  the  findings  of  Almen  and  Bruce  (Refs.  9  and  4  9)  and 
of  Brown  and  Davies  (Refs.  15  and  19),  who  noted  saturation  of  the  tar¬ 
get.  From  the  observations  made  by  these  investigators,  a  description 
of  the  phenomenon  of  saturation  may  be  formed.  In  the  presence  of 
sputtering,  the  target  collects  essentially  100  percent  of  the  ion  beam 
of  energy  above  approximately  3  kev  at  the  initial  application  of  the 
beam.  This  can  be  shown  by  using  Fig.  20  (Ref.  9)  together  with 
Fig.  21  (Refs.  9,  19,  and  50).  As  can  be  seen  in  these  figures,  the 
collection  efficiency  decreases  continuously  from  a  maximum  value  of 
unity  at  the  initiation  of  the  beam  collection  to  a  minimum  of  zero  at  a 
time  when  the  target  becomes  saturated  with  beam  material.  At  this 
time  any  new  molecule  striking  the  target  travels  to  some  finite  depth, 
depending  upon  its  initial  energy,  where  it  reaches  a  plane  of  very  high 
density.  Finding  no  empty  site  in  which  to  be  trapped,  it  is  either 
reflected,  subsequently  diffusing  back  from  the  target',  or  it  replaces 
an  already  trapped  molecule  which  diffuses  back  from  the  target.  Re¬ 
emission  of  the  "striking"  molecules  after  saturation  as  opposed  to  re¬ 
emission  of  already  trapped  molecules  is  indicated  for  noble  gases  in  the 
saturation  studies  of  Brown  and  Davies. 

The  observed  changes  in  range  distribution  of  collected  molecules  in 
the  noble -gas -bombarded  targets  suggests  that  self -sputtering  of  these 
molecules  does  not  occur  (Fig.  22);  they  are  simply  driven  deeper  into 
the  target  (Ref.  50);  that  is,  the  governing  mechanism  for  re-emission 
appears  to  be  collision  with  target  and  collected  ions  as  opposed  to  sput¬ 
tering  by  incoming  ions.  This  follows  from  the  previously  cited  observa¬ 
tion  on  re -emission. 

Unfortunately,  no  such  experimental  observation  has  been  made  for 
other  ions.  From  the  available  information,  a  slight  difference  is  evident 
for  other  ions  in  which  a  stronger  bond  occurs  between  the  ion  and  lattice, 
for  example,  Cs+  on  Cu.  It  is  believed  that  for  such  combinations  the 
density  distribution  of  collected  ions  in  the  saturated  target  would  exhibit 
a  higher  concentration  of  trapped  ions  near  the  target  surface*,  resulting 
in  a  combination  of  self-sputtering,  ion-target  sputtering,  and  ion  re- 
emission.  In  such  a  distribution,  a  high  concentration  of  collected  mole¬ 
cules  is  in  a  situation  that  is  favorable  to  self-sputtering  (Fig.  22).  This 
suggests  that  the  sputtering  ratio  for  a  given  ion-target  combination  changes 
continuously  until  the  target  is  saturated,  thereafter  becoming  constant  and 


*It  is  important  to  note  that  the  actual  range  distribution  for  all  mole¬ 
cules  has  the  same  form  even  in  the  saturated  target.  The  final  distribution 
of  trapped  molecules  in  the  saturated  target  is  the  property  which  is  variant 
in  form. 
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representative  of  the  ion-target  mixture  as  opposed  to  the  target-ion 
beam  combination  alone.  This  effect  has  been  noted  by  Almen  and 
Bruce. 


To  explain  the  previous  observations,  a  simplified  model  can  be 
developed.  For  a  given  target,  assume  that  a  finite  number  of  sites 
exists  wherein  an  ion  can  be  collected,  the  average  number  depending 
upon  the  depth  to  which  the  ion  can  penetrate.  Assuming  that  the  dis¬ 
tribution  of  these  sites  is  such  that  the  probability  of  the  occupation  of 
a  site  by  an  incoming  ion  is  dependent  only  upon  the  total  remaining 
number  of  sites  available  to  it,  the  collection  rate  of  beam  material  can 
be  expressed  by: 

k  <f)  —  (f)  l  ~  A  (29) 

Mi  I  (30) 


Dq  (31) 

where  D  is  dependent  upon  the  diffusion  capabilities  of  the  trapped 
molecule  and,  where  self- sputtering  occurs,  on  the  value  of  Ss . 


dq  r 

dt 

0  = 

Further,  assume  that 

A  = 


The  sticking  fraction  f  is  obtained  from: 

,  -  c(i  -  J-) 

where  qs  is  the  saturation  value  of  target  (gm). 
Combining  Eqs.  (29),  (31),  and  (32) 

IT  +  (D  +  q  " 

and  therefore 

q  = 

leading  to 

q  = 


<f>i 


D  + 


cjy  c 
<ls 


,  -(»♦£)' 
1  —  e 


+  -a- 


<Is 


(f)c 


-<f>c  (^-+-L- 
_  „  V  00  9 s. 


(32) 


(33) 


(34) 


(3  5) 


Equation  (3  5)  may  be  used  to  explain  most  of  the  saturation  phenom¬ 
ena  reported  in  the  experiments  of  Almen  and  Bruce  as  shown  in  Fig.  20. 
This  treatment  is  similar  to  that  of  adsorption  phenomena.  The  important 
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difference  can  be  recognized,  however,  by  comparing  the  magnitude  of 
the  sticking  probability,  accommodation  coefficient,  and  saturation 
value.  The  first  two  are  usually  smaller  for  surface  phenomena,  where¬ 
as  the  third  is  invariably  smaller  by  a  factor  of  2  to  1,  000  or  more, 
depending,  of  course,  upon  the  degree  of  penetration  of  the  target  by  the 
beam  particles.  Since  the  range  of  an  ion  is  a  linear  function  of  energy 
over  a  wide  range  of  application  (Refs.  9,  19,  50,  and  51),  the  satura¬ 
tion  value  can  be  expected  to  be  a  linear  function  of  energy  in  view  of 
the  niodel  adopted;  that  is,  the  number  of  sites  available  to  an  array  of 
impinging  ions  is  proportional  to  their  mean  depth  of  penetration.  This 
phenomenon  is  also  reported  by  Almen  and  Bruce.  In  comparing  Eq.  (3  5) 
with  the  experimental  data  available,  the  value  of  c  was  found  to  vary 
only  slightly  from  unity.  Saturation  values  for  various  targets  of  interest 
for  application  to  electrical  propulsion  test  facilities  are  plotted  in 
Fig.  23.  Almen  and  Bruce  have  shown  that  A,  the  re-emission  of  col¬ 
lected  atoms,  in  Eq.  (29)  is  effectively  zero  for  most  ion  target  combina¬ 
tions  at  room  temperature  when  the  energy  of  bombardment  is  above 
approximately  2  kev.  By  using  radioactive  ions  as  a  tracer,  A  was  meas¬ 
ured  by  setting  9$  equal  to  zero  after  a  given  bombardment  and  checking 
specific  radioactivity  of  the  target  as  a  function  of  time.  Observations 
lasting  up  to  10  days  showed  no  substantial  decrease  in  specific  radio¬ 
activity  for  a  large  number  of  combinations  which  usually  included  the 
heavier  ions  as  beam  material.  Such  targets  are  considered  permanent 
in  nature.  As  would  be  expected,  the  re-emission  rate  and  therefore  the 
diffusion  coefficient  D  increase  with  temperature.  Various  constants  and 
values  applicable  to  Eq.  (34)  have  been  tabulated  in  Table  2.  By  use  of 
these  constants,  certain  useful  calculations  are  possible.  For  example, 
the  time  required  for  a  target  to  reach  98  percent  of  its  saturation  value 


is  equal  to  four  time  constants, 4 


of  Eq.  (34),  which  is  attained 


almost  immediately  after  starting  the  ion  engine. 


4.3  BEAM  COLLECTION 

Below  a  certain  critical  energy,  it  appears  likely  that  an  unlimited 
quantity  of  any  solid  element  might  be  retained  on  the  surface  of  a  satu¬ 
rated  target.  This  phenomenon,  covered  in  some  detail  by  Almen  and 
Bruce,  occurs  only  in  the  presence  of  what  they  called  "zero  sputtering".  * 
They  proposed  three  basic  requirements  for  "zero  sputtering":  (1)  the 
self-sputtering  ratio,  that  is,  the  sputtering  ratio  of  ions  on  the  same 


*This  term  is  meant  to  imply  zero  target  material  sputtering.  Sput¬ 
tering  of  the  collected  material  may  continue  to  occur  with  some  sputter¬ 
ing  ratio  less  than  unity. 

\ 
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target  material,  must  be  less  than  unity,  (2)  the  initial  sputtering  of  the 
target  must  not  be  too  large,  and  (3)  the  target  temperature  must  be 
sufficiently  low  to  avoid  evaporation  of  the  deposit.  In  view  of  the  above 
treatment,  one  basic  requirement  which  comprises  all  of  the  above  three 
can  be  made.  The  collection  efficiency  for  the  beam  must  be  greater 
than  its  self- sputtering  ratio.  Where  this  is  true,  the  situation  then  pro¬ 
ceeds  as  follows:  Initially,  sputtering  occurs  at  the  target  in  a  normal 
way;  however,  the  target  becomes  saturated  and  a  deposit  of  the  ion 
species  forms  on  the  surface  of  the  target.  The  target  will  then  gain  in 
weight  according  to  the  following: 

- 0  (k  -  SB)  <36) 

If  k  >  Ss,  the  saturation  of  the  target  seemingly  never  occurs,  and  since 
the  deposit  eventually  grows  out  from  the  target  face,  the  original  target 
has  an  infinite  lifetime  and  remains  undamaged  after  its  initial  exposure 
to  the  beam.  (A  finite  time  is  required  for  the  target  saturation  and  sur¬ 
face  deposition  by  material  of  the  ion  species.  ) 

Self- sputtering  measurements  are  available  in  a  number  of  refer¬ 
ences,  but  specific  cases  of  interest  are  not  available  at  this  time.  The 
use  of  the  zero  sputtering  concept  in  the  resolution  of  problems  depends 
upon  the  ability  of  propellant-target  combination  to  support  the  proposed 
phenomena.  This  in  turn  will  require  some  form  of  beam  decelerating 
system  to  reduce  the  ion  energy  to  the  necessary  magnitude  such  that 

k  >  Sg. 

In  the  event  that  "zero  sputtering"  does  not  occur,  saturation  of  a 
target  results.  When  saturation  is  obtained,  continual  target  impingement 
results  in  self- sputtering,  diffusion,  or  a  mixture  of  the  two  as  a  mechan¬ 
ism  for  re-emission  of  molecules  of  the  ion  species  in  addition  to  sputter¬ 
ing  of  the  target  species. 


4.4  THE  TARGET  AS  A  SOURCE  OF  PARTICLES 

Based  on  the  results  of  the  previous  discussion,  if  k  <  Ss,  an  analyt¬ 
ical  treatment  for  the  accommodation  of  an  ion  beam  except  for  very  short 
tests  will  deal  only  with  the  phenomena  occurring  after  target  saturation. 
Where  k  <  Ss,  the  target  now  acts  as  a  source  of  mass  given  by 

Q  =  -L  (Mi  +  SMt)  <37> 


If  k  >  Ss,  the  source  after  initial  target  saturation  and  deposition  of 
a  few  layers  of  beam  material  will  be  given  by 

Q '  =  </>  [ss  +  (1  -  k)]  (38) 
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The  former  case  will  be  considered  first,  although  the  method  devel¬ 
oped  is  directly  applicable  to  either  case.  The  importance  of  the  source 
given  by  Eq.  (3  7)  depends  upon  the  size  of  the  ion  engine  and  duration  of 
the  test.  Figure  24  is  a  plot  of  Eq.  (3  7)  for  a  number  of  values  of  antici¬ 
pated  engine  test  conditions.  The  source  given  by  Eq.  (3  7)  must  be  com¬ 
pletely  accommodated  for  the  continuous  operation  of  a  facility  to  prevent 
contamination  of  the  test  cell. 


Of  particular  importance,  for  example,  is  the  flux  of  particles  in 
the  direction  of  the  engine  thruster  or  chamber  instrumentation.  Refer¬ 
ence  to  Fig.  25  gives  the  geometry  for  the  following  expression  for  the 
number  sputtered  per  incident  ion  in  the  direction  of  the  engine: 

S  (X)  =  So  [sA  (x)]  [sR  (X)]  (39) 

Available  data  indicate  that  in  the  energy  range  from  a  few  kev  to  ten  kev 
and  for  values  of  X  up  to  70  deg 

SA  (X)  =  cos  (X)  (40) 


where  n  >  0. 


Sr  (X) 


1 

cos  n  (X) 


(41) 


Thus, 


S  (X) 


s0 

(cos  X)"-1 


(42) 


The  following  conclusions  can  now  be  drawn.  If  the  product  SaSr  >  1, 
a  tilted  target  cannot  be  justified.  This  is  the  case  up  to  at  least  ten  kev. 
But  as  mentioned  previously,  at  higher  energies  Sj\  becomes  more  concen¬ 
trated  in  the  normal  direction.  Therefore,  if  the  effect  is  found  great 
enough,  the  majority  of  the  flux  can  be  aimed  at  a  suitable  spot  in  the  col¬ 
lector.  However,  at  the  present  time  no  such  data  are  available  to  justify 
this  technique.  It  should  be  noted  that  certain  exceptions  to  Eqs.  (40)  and 
(41)  exist,  but  all  exceptions,  such  as  platinum,  reported  in  the  literature 
at  this  time  are  impractical  for  use  in  testing. 


For  angles  X  less  than  70  deg  and  energies  above  a  few  kev  the  angu¬ 
lar  distribution  of  re-emitted  beam  particles  has  been  shown  to  be  cosine. 
Assuming  direct  line  of  flight  from  the  target  to  the  engine  thruster,  an 
approximation  of  flux  impingement  on  the  thruster  can  be  made  by  assum¬ 
ing  that  the  target  is  a  point  source  and  that  a  uniform  distribution  over 
the  engine  subtends  a  solid  angle  AO.  Neglecting  any  secondary  sputter¬ 
ing  and  all  but  direct  line  of  flight  particles,  the  flux  impingement  can 
be  expressed  as  (Fig.  25) 

P  =  SA  (X)SR  (X)  4^  +  —  -4^  cos  (*)  (43) 

e  2  77  e  2  n  ' 
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The  first  term  results  from  sputtered  particles,  whereas  the  second 
results  from  the  beam.  Finally,  the  critical  parameter  for  the  opera¬ 
tional  system  is  the  collection  efficiency  for  /3  on  vital  parts  of  the 
thruster  where  failure  might  result. 


4.5  PRESSURE  INCREASE  CAUSED  BY  ENGINE  EFFLUX 


To  this  point,  material  added  to  the  system  has  been  characterized 
only  by  the  number  of  molecules  and  their  mass.  For  certain  applica¬ 
tions  this  would  be  sufficient  and  necessary  parameters  could  be  deter¬ 
mined  experimentally.  However,  to  determine  the  system's  ability  to 
accommodate  such  material,  it  sometimes  becomes  necessary  to  con¬ 
sider  also  the  mean  velocity  of  particles  and  the  effect  of  reflection 
from  various  points  in  the  collector.  One  such  approach  is  to  determine 
the  equivalent  pressure  increase  in  a  chamber  caused  by  the  new  material 
which  emerges  from  the  target  upon  the  introduction  of  an  ion  beam  into 
the  system.  One  important  result  of  this  treatment  is  that  it  displays  the 
relative  importance  of  the  beam  and  the  sputtered  particles  on  test  cham¬ 
ber  conditions.  The  method  of  Mickelsen  and  Childs  was  modified  and 
extended  for  this  application.  To  summarize  their  method,  an  equilibrium 
situation  is  assumed  wherein  the  total  momentum  transfer  to  the  walls  of 
the  collector  is  given  as  a  summation  over  n  collisions,  n  being  the  num¬ 
ber  required  for  the  condensation  of  beam  material  on  a  cold  collector 
surface.  If  the  pressure  experienced  at  the  target  is  neglected,  the  pres¬ 
sure  contribution  within  the  test  chamber  by  particles  making  the  nth  col¬ 
lision  is  given  by 


APn  = 


Vp 

2A„  L 


1  +  (1  -  k)(l  -  a) 


x][(i  - 


k) (1  -  a) 


% 


(44) 


where  the  zeroth  collision  is  that  made  on  the  target.  Since  the  accom¬ 
modation  coefficient  and  collection  efficiency  will  in  most  cases  of  interest 
be  very  near  unity,  only  the  first  few  collisions  are  of  importance;  that  is, 
the  series  converges  very  rapidly.  For  this  reason,  it  is  easily  shown 
that  very  little  error  is  introduced  by  allowing  n  to  approach  infinity. 
Hence, 


A  Pb 


4> b  v. 

2  A 


1  +  (1  -  k )  ( 1 


-■>*]£ 

Jn  =  1 


(1  -  k)(l  -  a) 


(45) 


r 


=  ^bVo  l 
2  A 


[  1  +  (1  -  k)(  1  -  a)‘4  ][  (I  -  k  )  (  1  -  a/4  ] 
1  -  (  1  -  k  )  (  1  -  a  )  ‘4 


a  random  angular  emission  distribution  results  in  the  factor  of  one  half. 
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The  same  result  is  derived  if  n  collisions  are  considered  and 

r  i n 

(1  -  k)(l  -  a)V2J  »  0  (46) 

Pressure  on  the  target  face  (the  zeroth  collision)  being  neglected,  the 
summation  begins  with  unity. 


From  previous  discussions,  it  is  concluded  that  the  collision  with 
the  target  cannot  be  considered  similar  to  remaining  collisions  made 
with  the  walls  of  the  collector  for  energy  and  saturation  reasons.  Inter¬ 
actions  at  points  in  the  collector  system  other  than  at  the  target  face 
will  be  surface  interactions,  whereas  beam  target  interactions  are  not 
strictly  surface  interactions.  In  the  latter  case,  the  impinging  mole¬ 
cules  penetrate  through  a  number  of  lattice  layers  before  they  are  re- 
emitted  from  the  target;  therefore,  the  accommodation  coefficient  of 
these  beam  particles  is  expected  to  be  much  greater  at  the  target  than 
at  any  other  collision  surface,  and  the  average  velocity  of  beam  particles 
emerging  from  the  target  is  given  by 

V  =  Vo(l  -  at/4  (47) 

The  magnitude  of  V  was  discussed  previously  in  section  2.  6.  2. 


Modification  of  Eqs.  (44)  and  (45)  to  consider  the  target  itself  as  a 
source  of  particles  leads  to  a  different  expression  for  the  pressure  con¬ 
tribution  of  the  particles  after  leaving  the  target.  A  general  angular 
distribution  will  be  assumed  for  the  present  by  inserting  a  constant 
factor  ip  in  calculations.  For  the  nth  collision,  considering  first  the 
beam  particles, 


APb 


iAb0b  Vq  (i  -  *t)1/2 

A, 


[(1  -  k)(l  -  a ) 1/2 


n  “  1  r 


1  + 


(l  -  k)(l  -  a ) 1/2  ]  (48) 


A  summation  over  n  leads  to 


APb 


lph<t>b  Vo  (1  -  [  (1  -  k)(l  -  a)  *4  +  1  ] 

A0  |  [  1  -  (1  -  k  )  ( 1  -  a)l4  ] 


(49) 


Until  now  no  mention  has  been  made  of  the  second  term  of  the  right 
hand  side  of  Eq.  (37),  which  describes  the  source  of  the  sputtered 
particles.  Similarly,  the  previous  treatment  (Eq.  49)  for  primary  beam 
particles  can  be  applied  to  the  sputtered  particles;  however,  the  results 
must  be  properly  interpreted  for  each  case  before  conclusions  of  any  value 
can  be  made,  especially  for  certain  ion-target  combinations.  This  is  true 
because  the  collection  efficiency  for  some  of  the  combinations  to  be  con¬ 
sidered  is  near  unity.  In  this  case,  the  angular  distribution  of  emergent 
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particles  is  very  important  since  there  is  a  large  pressure  variation  in 
the  collector;  that  is,  the  flux  of  material  is  concentrated  upon  a  local¬ 
ized  portion  of  the  large  vacuum  chamber.  For  the  same  reason,  pres¬ 
sure  readings  in  such  a  chamber  must  be  properly  interpreted  since 
these  atoms  could  only  be  detected  by  a  properly  oriented  "nude"  sens¬ 
ing  element.  Such  atoms  would  plate  out  on  any  gage  tubulation  before 
entering  the  electrodes  and  the  resulting  gage  readings  would  be  erron¬ 
eously  low.  Where  sticking  fractions  are  lower,  this  problem  is  not 
encountered.  The  following  additional  assumptions  are  required  to  ex¬ 
tend  the  theory  to  consider  sputtered  particles  effects: 

1.  Sputtering  at  any  point  in  the  collector  system  other  than 
at  the  target  is  considered  negligible.  This  is  true  be¬ 
cause  a  very  high  accommodation  coefficient  for  a 
normally  energetic  beam  results  in  low  ion  energies  for 
the  re-emitted  particles.  In  the  case  of  normal  incident 
beams,  the  primary  sputtered  particles  are  so  low  (20  to 
50  ev)  that  they  can  induce  very  little  secondary  sputtering. 

The  same  cannot  be  said  for  an  obliquely  incident  beam. 

An  increase  in  reflection  of  high  energy  particles  with 
very  little  accommodation  can  result  in  secondary  sputter¬ 
ing.  Even  this  effect  will  be  shown  negligible,  however, 
for  angles  of  target  rotation  as  high  as  60  deg. 

2.  The  remainder  of  the  collector  is  not  saturable.  This 
assumption  follows  if  the  surface  temperature  of  the  col¬ 
lector  walls  is  such  that  condensation  can  occur.  The 
relatively  low  energies  of  the  particles  result  in  their  con¬ 
densation  and  accumulation. 

3.  Gettering  can  be  considered  separately.  The  possibility  has 
been  recognized  that  sputtering  phenomena  might  be  harnessed 
to  augment  the  pumping  components  of  the  test  chamber.  Yonts 
has  noted  the  observation  of  a  reduction  in  residual  gas  pres¬ 
sure  when  the  materials  are  such  that  gettering  is  possible 
(for  example,  Ca+  on  Cu).  For  materials  not  conducive  to 
gettering,  a  pressure  rise  was  noted  (Ref.  32).  In  any  event, 
the  following  treatment  is  independent  of  gettering  of  the 
residual  gases,  and  a  further  study  will  be  necessary  to  in¬ 
clude  this  phenomenon. 


Using  the  same  treatment  for  primary  beam  particles,  the  following 
expression  of  pressure  contributions  in  the  nth  collision  for  sputtered 
particles  is  obtained: 


AP„  = 


<ftb  Vs  S 


(1  -  k.)(l 


1/  -|  n  —  1  [- 

0/2J 


1  +  (1  -  k)(l  _  a) 
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A  summation  over  all  collisions  yields 

y  Ap  y>s <ftb  V,  s  [  1  +  (l  -  k')(i  -  a)^  ] 

n  =  l  S  Ac  [  1  -  (  1  -  k  )  (  1  -  a  )%  ] 

To  find  the  total  pressure  caused  by  the  operation  of  an  engine,  add 
Eqs.  (49)  and  (51). 


AP 

If  it  is  assumed  that 

and  thus 


£  ( A  Ps  +  APb) 

n  =  1 


k  =  a  «=  1 

(1  -  k)(l  -  a/4  »  0 


(52) 


(53) 

(54) 


assuming  a  cosine  distribution  for  the  total  efflux,  making  the  factor 
xjj  =  2/3,  the  pressure  in  the  chamber  may  be  estimated  by: 

AP  =  [SMSVS  +  MbVo(l  ~  at ) ^ 


(55) 


Equation  (55)  is  plotted  in  Fig.  26  and  compared  with  experiments 
conducted  by  Richley  and  Cybulski.  While  Eq.  (55)  is  derived  from  a 
highly  idealized  model,  it  indicates  what  parameters  are  involved  in  a 
determination  of  pressure  changes  and  how  such  pressure  changes  vary. 


The  relative  importance  of  sputtered  and  re-emitted  particles  from 
the  target  may  be  found  by  comparing  their  contributions  to  the  change 
in  chamber  pressure  after  target  saturation  occurs.  This  may  be  ob¬ 
tained  by  taking  the  ratio  of  Eq.  (51)  with  Eq.  (49).  The  sputtering  ratio 
is  obtained  from  Eq.  (28),  and  it  is  assumed  that  Eq.  (54)  is  valid  for 
both  beam  and  sputtered  target  materials.  We  obtained  the  following 
relationship: 

A  Ps  =  Smax  (1  ~e~~CEi)  V  Ms  e5 
APb  VMi(l~at)Ei 

The  greatest  uncertainty  lies  in  the  proper  value  for  (1  -  a-^Ei. 

The  upper  and  lower  limits  can  be  established  based  upon  the  follow¬ 
ing  argument.  Since  the  mechanism  for  re -emission  of  beam  particles 
is  essentially  the  same  as  that  for  sputtered  particles,  the  upper  limit 
for  their  energy  is  conceivably  that  of  the  sputtered  particles.  However, 
since  these  beam  particles  penetrate  the  target  and  are  thermalized  to 
some  degree,  their  lower  energy  limit  can  be  estimated  as  that  energy 
corresponding  to  the  localized  target  temperature.  Until  further  studies 
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are  conducted  in  this  area  the  following  limits  are  assumed: 

-f-  KTt  <  (1  -  at)Ei  <  Es  (E 

Equation  (56)  is  plotted  (Fig.  2  7)  for  typical  values  of  the  limits  in 
Eq.  (57).  It  can  be  seen  from  Fig.  27  that  the  effects  of  sputtered  par¬ 
ticles  are  much  greater  than  those  caused  by  the  re -emitted  beam 
particles. 


4.6  COLLECTOR  SYSTEM  PROBLEMS 

The  effects  of  the  various  phenomena  inherent  in  the  collection  and 
termination  of  an  energetic  beam  upon  the  operational  capability  of  an 
electrical  propulsion  test  facility  will  now  be  considered.  For  long 
term  tests,  a  major  problem  is  the  lifetime  of  the  collector  system. 

The  rate  of  erosion  of  a  collector  surface  may  be  described  by  the  fol¬ 
lowing  expression  (assuming  a  uniform  current  distribution  impinging 
normally  upon  the  target): 

rate  of  erosion  =  — ^ —  (58) 

Ne 


Figure  28  presents  a  series  of  curves  of  the  erosion  rate  per  incid¬ 
ent  current  density  for  mercury  and  cesium  ions  on  aluminum,  iron, 
copper,  and  nickel  targets.  It  is  clear  that  for  long  term  engine  tests, 
the  conventional  metallic  plate  collector  may  not  be  used  unless  modi¬ 
fications  are  made.  For  example,  a  mercury  ion  beam  having  a  current 
density  of  30  ma/cm.2  with  a  beam  energy  of  6  kv  will  erode  approximately 
2  cm  of  target  material  after  30  hr  of  bombardment.  Electrical  propul¬ 
sion  systems  operating  for  extended  periods  of  time  (1,  000  hrs  or  longer) 
could  erode  meters  of  target  material.  Considering  the  nature  of  applica¬ 
tion  of  electrical  propulsion  systems,  the  necessity  of  tests  of  such  length 
is  evident  (Ref.  1),  and  the  sputtering  problem  demands  resolution  before 
long-term  test  programs  can  be  performed  in  ground  test  facilities. 

For  tests  of  a  shorter  duration,  the  erosion  rate  will  not  be  an  import¬ 
ant  factor.  However,  the  particle  efflux  defined  by  Eq.  (37)  introduced 
into  the  chamber  because  of  the  operation  of  the  engine  will  limit  the  degree 
of  simulation  possible  by  the  ground  test  facility.  In  a  two-fold  problem 
for  long  term  tests,  the  lifetime  of  the  system  as  well  as  the  efflux  prob¬ 
lem  must  be  resolved  before  an  effective  collector  system  is  achieved. 

4.6.1  Solid  Collectors 


The  solid  target  collector  is  restricted  to  low  total  ion  flux  tests 
except  under  one  condition  to  be  mentioned  presently.  Therefore,  the  use 
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of  Fig.  2  8  lends  a  rough  estimate  of  the  capabilities  of  such  a  system 
where  the  unmodified  beam  impinges  normally  on  such  a  target.  Exam¬ 
ination  of  the  relations  covered  in  past  sections  indicates  that  the  most 
favorable  conditions  for  the  collector  are  obtained  at  low  energies  (up 
to  a  few  kev)  for  a  normally  incident  beam. 

A  search  for  optimum  collector  geometry  has  been  conducted 
(Ref.  46).  Some  of  the  configurations  include  honeycomb-type  surfaces 
and  variously  arranged  fin- covered  walls.  Preliminary  data  from  such 
arrangements  indicate  that  very  little  is  gained  by  the  introduction  of 
fins  and  wells  such  as  an  egg  crate  structure  into  the  collector.  This 
can  be  understood  by  examining  Eq.  (55).  The  value  for  area  to  be  used 
in  that  equation  must  be  an  unshadowed  area  and  hence  can  be  considered 
to  be  the  area  of  an  imaginary  smooth  surface  surrounding  the  target. 

The  magnitude  of  such  a  surface  cannot  be  increased  by  the  addition  of 
fins  or  plates  because  the  first  collision  with  the  wall  contributes  by  far 
the  most  to  the  pressure  (k  =  1).  Such  fins  can  only  serve  as  cryosur- 
faces  to  reduce  residual  gas  pressure  or  for  beams  composed  of  mole¬ 
cules  which  do  not  condense  in  the  first  few  collisons.  Furthermore,  the 
use  of  finned  or  honeycombed  surface  targets  does  not  prove  feasible  for 
long  term  tests  where  the  beam  energy  is  not  reduced  because  of  surface 
erosion  problems. 

4.6.2  Long  Term  Tests  Using  Solid  Collectors 


For  long  term  tests,  Fig.  28  illustrates  that  solid  targets  could  be 
used  in  an  unmodified  beam  only  by  continuously  replacing  targets  as 
they  are  deteriorated  by  the  beam.  Replacing  a  large  target  without 
breaking  the  vacuum  is  quite  a  formidable  task. 

One  very  important  phenomenon  which  has  been  observed  may  offer 
the  solution  to  problems  presented  by  sputtering.  It  was  previously 
described  as  "zero  sputtering"  and  can  be  controlled  and  employed  only 
if  the  beam  can  be  decelerated  by  some  means  to  an  energy  level  so  that 
the  collection  efficiency  at  the  target  is  greater  than  the  self-sputtering 
coefficient.  Then  equilibrium  conditions  would  be  established  without 
appreciable  damage  to  the  target  for  tests  of  any  duration. 

In  principle,  it  is  possible  to  retard  the  beam  by  having  the  collector 
at  a  high  voltage.  However,  the  retarding  field  also  acts  as  an  electro¬ 
static  lens,  which  causes  defocusing  of  the  beam  at  low  energies.  By 
using  a  suitably  shaped  collector  arrangement,  Elbeck,  Oleson,  and 
others  have  reported  a  solution  to  overcome  this  defocusing  (Ref.  54). 

For  any  decelerating  system,  a  convenient  method  for  ascertaining 
the  desired  ion  energy  is  similar  to  that  used  by  Wehner  to  determine 


27 


AEDC-TDR-64-105 


sputtering  thresholds  (Ref.  55)  employing  emission  spectroscopy.  When 
zero  sputtering  conditions  are  met  and  equilibrium  is  established,  no 
spectral  lines  characteristic  of  the  target  can  be  detected  by  the  spec¬ 
trometer.  If  it  is  possible  to  achieve  the  phenomenon  described  as  zero 
sputtering,  there  is  little  question  that  herein  lies  the  solution  of  the 
sputtering  problem  in  electrical  propulsion  test  facilities. 

Some  of  the  criteria  which  must  be  considered  in  the  selection  of 
materials  for  collector  systems  are: 

1.  electrical  conductivity 

2.  thermal  conductivity 

3.  sputtering  coefficient 

4.  secondary  emission  properties 

5.  vapor  pressure  at  operating  temperature  of  thruster 
components 

6.  economy 

The  above  criteria  and  the  data  in  Figs.  29a,  b,  c,  and  d  (Refs.  56 
and  57)  on  the  sputtering  yield  of  various  targets  by  Hg+  and  Xe+  indicate 
that  the  most  favorable  target  materials  are  aluminum,  nickel,  and  iron. 
The  Xe+  ion  is  chosen  because  it  behaves  quite  similarly  to  Cs+  (Ref.  20). 
Even  though  these  metals  do  not  exhibit  the  lowest  sputtering  yields,  they 
were  chosen  for  economy  or,  more  important,  for  the  high  vapor  pres¬ 
sure  characteristics  they  exhibit  at  the  operating  temperatures  of  typical 
thruster  components.  For  example,  the  components  of  cesium  contact 
engines  generally  operate  at  temperatures  of  approximately  1400°K. 
Therefore,  while  a  material  such  as  carbon  appears  because  of  its  small 
sputtering  ratio  to  be  favorable,  its  use  must  be  avoided  because  its  vapor 
pressure  at  1400°K  is  very  low  (approximately  10"  10  torr)  and  could  cause 
failure  by  condensing  on  vital  parts  of  the  ion  engine  (for  example  on  the 
ionizer).  Of  the  three  collector  materials  mentioned  above,  aluminum 
possesses  the  most  favorable  overall  properties  (note  the  erosion  rates 
in  Fig.  28). 

4.6.3  Liquid  Collector  Systems 


The  use  of  a  liquid  metal  as  a  collector  target  has  been  considered 
(Ref.  1)  and  appears  to  offer  some  advantages  over  the  conventional  solid 
system.  It  may,  for  example,  eliminate  collector  system  lifetime  prob¬ 
lems.  In  addition,  some  liquid  metals  have  desirable  characteristics  when 
used  as  gettering  materials.  Furthermore,  this  suggests  the  possibility 
of  collecting  and  reusing  expended  propellant  by  using  the  same  type  of 
material  for  the  liquid  metal  target  as  is  used  for  the  propellant. 
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The  practicality  of  a  liquid  metal  system  is  being  investigated  and 
comparisons  made  with  present  day  solid  systems.  Experiments  are 
required  to  determine  if  a  significant  difference  exists  in  the  sputtering 
of  a  liquid  and  solid  target.  Some  work  has  been  reported  at  ion  ener¬ 
gies  below  one  kev  (Ref.  57)  (Figs.  30  and  31);  however,  no  precise 
statements  can  be  made  about  sputtering  ratios,  distribution,  or  ener¬ 
gies  of  sputtered  particles  for  higher  ion  energies. 


5.0  SUMMARY  AND  CONCLUSIONS 


The  foregoing  study  comprises  a  compilation  and  investigation  of 
studies  on  sputtering  made  by  numerous  investigators  to  date.  This  was 
done  with  the  objective  of  obtaining  information  necessary  for  the  design 
of  space  simulation  chambers  for  testing  advanced  propulsion  devices 
for  which  exhaust  beam  termination  represents  a  serious  problem  in 
view  of  its  ablative  action  upon  the  target  and  chamber  walls  and  its  pos¬ 
sible  interaction  with  the  system  to  be  tested.  Since  no  analytical  treat¬ 
ment  was  found  in  which  sputtering  and  saturation  phenomena  were  accu¬ 
rately  described  over  a  broad  energy  range,  a  theoretical  study  was  made. 

Equations  were  derived  leading  to  the  sputtering  ratio  characteristic 
of  two  types  of  collision  models  which  are  applicable  in  different  energy 
ranges.  The  resulting  equations  then  serve  as  the  upper  and  lower 
boundaries  for  a  general  expression  throughout  the  two  energy  ranges 
applicable  to  electrostatic  propulsion.  The  form  of  this  equation  com¬ 
pares  favorably  with  experimental  results  obtained  at  other  laboratories. 

Also  an  attempt  was  made  to  establish  a  theoretical  basis  for  satura¬ 
tion  phenomena.  The  resulting  equations  are  in  good  agreement  with 
experimental  data.  Based  on  the  above  study  the  following  conclusions 
are  drawn: 

1.  To  describe  the  sputtering  phenomenon  in  the  region  of  interest 
in  ion  propulsion,  at  least  two  different  collision  models  are  required. 

2.  The  termination  of  an  ion  beam  at  a  target  can  result  in  the 
addition  of  free  target  and  ion  beam  particles  into  the  test  chamber. 

The  amount  of  contamination  is  shown  to  be  appreciable  if  a  collector 
is  not  designed  to  accommodate  this  material. 

3.  Positive  and  negative  ion  emission  from  the  collector  is  much 
less  than  secondary  electron  emission,  which  in  turn  is  much  less  than 
neutral  particle  emission. 

4.  Under  operating  conditions,  a  solid  target  will  become  saturated 
with  beam  particles  within  a  few  milliseconds  except  in  the  presence  of 
"zero  sputtering". 
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5.  If  the  energy  of  the  ion  beam  can  be  controlled  without  altering 
operating  conditions,  it  may  be  possible  to  collect  more  material  than 
is  sputtered,  thus  eliminating  target  lifetime  problems. 

6.  It  appears  that  xenon  sputtering  data  can  be  substituted  where 
cesium  data  are  not  available. 

7.  Under  long  term  test  conditions  where  solid  targets  are  not 
adequate,  the  use  of  liquid  metal  targets  appears  feasible. 

8.  Experimental  work  is  necessary  in  the  following  areas; 

a.  Sputtering  and  saturation  of  solid  and  liquid 
targets  by  prospective  propellant  materials. 

b.  The  effect  upon  flight  simulation  of  beam  decel¬ 
eration,  sputtered  particle  flux,  and  secondary 
emission. 
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TABLE  1 

J3 

VALUES  OF  Smax 

AND  FROM  EQ 

'.  (28) 

Target 

K 

t  ■( 

■'tt-a 

Ion 

Material 

Snaax 

(kev)_l 

ev 

Hg+ 

Cu  (Ref.  58) 

14 1 0 

0.  179 

60 

Hg+ 

Ni  (Ref.  5  8)  V 

8.  0 

0.  218 

80 

Hg+ 

Fe  (Ref.  58) 

6.  0 

0.  179 

65 

Hg+ 

A1  (Ref.  45)  t 

2.4 

0.  218 

130 

Cs  + 

Cu  (Ref.  32)  Aj 

15 

0.  34 

60 

TABLE  2 

SATURATION  VALUES  AND  TIME  CONSTANTS 


Hgm  cm-2 

cm2  sec-l  amp-*  x  10^ 

Ion 

Target 

qs  (45  kev) 

Calculated  decay  constant 
current  density 

Kr+ 

A1 

5.2 

1.  6 

Fe 

2.  1 

4.0 

Ni 

1.  6 

5.  2 

Xe  + 

A1 

4.2 

3.  1 

Fe 

1.4 

9.  5 

Ni 

1.  5 

8.  8 
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Fig.  1  Typical  Sputtering  Yield  Curve 
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Fig.  2  Nickel  Sputtering  Yields  at  150  ev  vs  Argon  Pressure 
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Normal  Incidence 

b.  Fe 

Fig.  4  Continued 
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Fig.  4  Continued 
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d.  Cu,  W,  Mo 
Fig.  4  Continued 
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Fig.  4  Continued 
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Krypton  Ions,  45  kev 


f.  Ag,  Cu,  Ni,  To 
Fig.  4  Concluded 


45 


03 


M  +  .  ocn  Relative  Intensity 
Hg  at  250  ev  J 


Fig.  5  Angular  Distribution  of  Sputtered  Particles  at  Oblique  Incidence  (250 ev) 
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Fig.  7  Angular  Distribution  of  Sputtered  Particles  at  Oblique  Incidence  (20  kev) 
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a.  Hg+  on  Ni 

Fig.  8  Angular  Distribution  of  Sputtered  Particles 
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b.  Hg+  on  Fe 
Fig. 8  Continued 
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Fig.  9  Angular  Distribution  of  Sputtered  Particles 
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Fig.  10  Velocity  Distributions  of  Sputtered  Copper  Atoms 


AEDC-TDR-64-105 


CJI 


7 


6 


5 


Lf\ 

I 

o 


x  4 

O 


o 


2 


1 


0 


Sym 

□ 

A 

o 

ncidence  Angle,  deg 

W-60 

W-  0 

Mi-  n 

X 

Hg+  Inch 

Cu- 

jent  Ion 

0 

Refs.  17 

and  18 

('  " 

0  100  200  300  400 


Fig.  11  Average 


SB] 


SHIil 


SB 


Ion  Energy,  ev 

ity  of  Sputtered  Particles  v 


s  Ion  Energy 


1 


R  -  64 - 1 05 


0  0.4  0.8  1.2  1.6  2.0  2.4  2.8 

Energy,  kev 


Fig.  14  Effect  of  Temperature  and  Surface  Condition  on  Secondary  Emission  (Cs+) 
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Fig.  15  Sputtering  Yields  for  Different  Ions  ot  45  kev 
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Fig.  18  Sputtering  Yield  vs  Ion  Energy  (Mercury  Bombarding  Copper  and  Nickel) 
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Fig.  19  Sputtering  Yield  vs  Ion  Energy  (N2+  on  Ni) 
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Fig.  24  Source  of  Material  from  Saturated  Target  vs  Sputtering  Ratio 
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Fig.  27  Ratio  of  Sputtered  Particle  Pressure  to  Beam  Pressure  vs  Ion  Energy 
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Fig.  31  Sputtering  Yields  for  Liquid  Tin  and  Gallium  (Hg+) 
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